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Abstract 

We present the construction of an exponentially accurate time-dependent Born- 
Oppenheimer approximation for molecular quantum mechanics. 

We study molecular systems whose electron masses are held fixed and whose nuclear 
masses are proportional to e~ 4 , where e is a small expansion parameter. By optimal 
truncation of an asymptotic expansion, we construct approximate solutions to the 
time-dependent Schrodinger equation that agree with exact normalized solutions up 
to errors whose norms are bounded by C exp ( — 7/e 2 ) , for some C and 7 > 0. 
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1 Introduction 



In this paper we construct exponentially accurate approximate solutions to the time-depen- 
dent Schrodinger equation for a molecular system. The small parameter that governs the 
approximation is the usual Born-Oppenheimer expansion parameter e, where e 4 is the ratio of 
the electron mass divided by the mean nuclear mass. The approximate solutions we construct 
agree with exact solutions up to errors whose norms are bounded by C exp ( —7 /e 2 ) , for 
some C and 7 > 0, under analyticity assumptions on the electron Hamiltonian. 

The Hamiltonian for a molecular system with K nuclei and N — K electrons moving in 
I dimensions has the form 

K 4 N -1 

j=l * m j j=K+l £U H i<j 

Here Xj G IR l denotes the position of the j th particle, the mass of the j th nucleus is £~ 4 Mj 
for 1 < j < K, the mass of the j th electron is rrij for K + 1 < j < N, and the potential 
between particles i and j is Vij. For convenience, we assume each Mj = 1. We set d = Kl 
and let X = (Xi, X 2 , . . . , Xk) G IR d denote the nuclear configuration vector. We can then 
decompose H(e) as 

H(e) = -yA x + h(X). 

The first term on the right hand side represents the nuclear kinetic energy, and the second 
is the "electron Hamiltonian" that depends parametrically on X. For each fixed X, h(X) is 
a self-adjoint operator on the Hilbert space H e \ = L 2 (M^ N ~ K ^ 1 ). 

The time-dependent Schrodinger equation we approximately solve in L 2 (IR d , 7i e \) as 
e -> is 

ie2 ^t = ~^ A x^ + h(X)^ 

Asymptotic expansions in powers of e of certain solutions to this equation are derived in 
0, D, p. We obtain our construction by truncating these expansions after an e-dependent 
number of terms, in an effort to minimize the norm of the error. Similar strategies have been 
used to obtain exponentially accurate results for adiabatic approximations BR [TSL [19], |16[ 



and semiclassical approximations |L4], |To] , both of which play roles in the Born-Oppenheimer 
approximation we are studying here. 

Roughly speaking, the time-dependent Born-Oppenheimer approximation says the fol- 
lowing for small e: The electrons move very rapidly and adjust their state adiabatically as the 
more slowly moving nuclei change their positions. If the electrons start in a discrete energy 
level of h(X), they will remain in that level. In the process, the electron states create an 
effective potential in which the motion of the heavy nuclei is well described by a semiclassical 
approximation. The asymptotic expansions show that this intuition is valid up to errors of 
order e k for any k. 

In Born-Oppenheimer approximations, adiabatic and semiclassical limits are being taken 
simultaneously, and they are coupled. Analysis of errors for the adiabatic and semiclassical 
approximations shows that they are each accurate up to errors whose bounds have the form 
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C exp( — j/e ) (I7|, [14]]. Non-adiabatic transitions are known in some systems to be of 
this order, and tunnelling in semiclassical approximations makes contributions of this order. 
Thus, one cannot expect to do better than approximations of this type because of two 
well-known physical phenomena that Born-Oppenheimer approximations do not take into 
account. 

In some systems, tunnelling might dominate the error. In some, non-adiabatic electronic 
transitions may dominate. In others, the two effects can be of comparable magnitude. 

One of the motivations for our work is to generate a "good" basis upon which to build a 
"surface hopping model" that would accurately describe non-adiabatic electronic transitions. 
Prior authors (see, e.g., [|^, |28], [29|, ||) have proposed such models based on the zeroth order 
time-dependent Born-Oppenheimer approximation. Using the zeroth order states as a basis 
of the surface hopping model, the non-adiabatic transitions appear at order e 2 . This is 
huge compared to the exponentially small physical phenomenon one would like to study, 
and we believe interference between transitions that occur at different times is responsible 
for the exponential smallness of the physically interesting quantity. Our view is that by 
choosing a much better set of states on which to base the model, one will obtain a much 
more useful approximation. Sir Michael Berry |3], f|, has advocated such ideas for the 
somewhat simpler adiabatic approximation (which does not have the complications of the 
nuclear motion). These ideas have been used in [L8| to prove the accuracy of certain 



results for non-adiabatic transitions that are exponentially small. 
Remarks: 

1. There are some other exponentially accurate results in the general topic of Born- 
Oppenheimer approximations. The prior results come from study of the time-independent 
Schrodinger equation and depend on global properties of the system. Our results are time- 
dependent and make use of local information. 



Klein |2(J and Martinez [^, |23|, [24| show that resonances associated with predissociation 
processes have exponentially long lifetimes. Benchaou and Martinez |1|, |2[] also show that 
certain S-matrix elements associated with non-adiabatic transitions are exponentially small. 

2. The papers cited in the previous remark obtain estimates that depend on the global 
structure of the electron energy levels. The results we obtain depend on a particular classical 
path. When the path stays away from the nuclear configurations where the gap between 
relevant electonic levels is minimized, one would expect the non-adiabatic errors from our 
approximation to be smaller, i.e., both results would obtain errors of order exp(— r/e 2 ), but 
we would obtain a larger value of T. 

We expect this because in our case, the Landau-Zener formula predicts that our T should 
come from the minimum gap between eigenvalues on the classical path, rather than the global 
minimum gap. 

3. From a mathematical point of view, the optimal truncation procedure in this context was 
first stated for the adiabatic approximation for two component systems of ODE's by Berry 
0, |||. It was first proved to yield exponentially accurate results for Hilbert space valued 
ODE's by Nenciu p5| . In [|14], |15|] we used this idea for the semiclassical approximation, 
which is a complex valued PDE setting. The present paper can be viewed as extending these 
ideas to a Hilbert space valued PDE setting. 
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1.1 Hypotheses 

We assume that the electron Hamiltonian h(X) satisfies the following analyticity hypotheses: 

Ho (i) For any X G lR d , h(X) is a self-adjoint operator on some dense domain D C Tt e i, 
where 7i e \ is the electronic Hilbert space. We assume the domain T> is independent 
of X and h(X) is bounded from below uniformly in IR d . 

(ii) There exists a 5 > 0, such that for every ip G T>, the vector h(X)ip is analytic in 
S s = {ze(U d : \1m( Zj )\<6, j = l,...,d}. 

Hi There exists an open set H C lR d , such that for all X G H, there exists an isolated, 
multiplicity one eigenvalue E(X) of associated with a normalized eigenvector 

G "H e i- We assume without loss that the origin belongs to H. 



Remarks: 1. Hypothesis Ho implies that the family of operators {h(X)}x£S s is a holo- 
morphic family of type A. 

2. It follows from H and H x that there exists 5' G (0, 5) and ? C H such that the 
complex and vector valued functions E(-) and $(•) admit analytic continuations on the set 
E 5 , = {ze(U d : Re(z) G S' and \lm( Zj )\ < 5', j = 1, . . . , d}. 

1.2 Summary of the Main Results 



Our main results are stated precisely as Theorem [4.1| in Section ^. Two generalizations of 
this result are presented in Section 0. 



Roughly speaking, Theorem ^4.1| states the following: 



Under hypotheses H and Hi, we construct ^*(X, t, e) (that depends on a parameter g) 
for t G [0, T]. For small values of g, there exist C(^) and r(g) > 0, such that in the limit 



e-^)/^(X,0,e) -*,(X,t,c) r2/iDd ^ N < C(<?) e 



i 2 (R d ,H el ) 



In the state t,e), the electrons have a high probability of being in the electron state 

$(X). For any b > and sufficiently small values of g, the nuclei are localized near a classical 
path a(t) in the sense that there exist c(g) and 7(g) > 0, such that in the limit e — > 0, 



1/2 

7(9)/e 2 



i(t)|>6 

The mechanics of the nuclear configuration a(t) is determined by classical dynamics in the 
effective potential E(X). 

Two theorems in Section |8] generalize this result. The first allows the time interval to 
grow as e tends to zero. The second allows more general initial conditions. 
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2 Coherent States and Classical Dynamics 



In the construction of our approximation to the solution of the molecular Schrddinger equa- 
tion, we need wave packets that describe the semiclassical dynamics of the heavy nuclei. In 
the present context, the semiclassical parameter is h = e 2 . We make use of a convenient set 
of coherent states (also called generalized squeezed states), that we express here in terms of 
the semiclassical parameter h. 

We recall the definition of the coherent states (pj(A, B, h, a, rj, X) that are described in 
detail in |I3fl . A more explicit, but more complicated definition is given in fT2|j . 

We adopt the standard multi-index notation. A multi-index j = (ji, j 2 , . . . , jd) is a 
rf-tuple of non-negative integers. We define \j\ = J2f=i 3k, X^ — X^X^ 2 ■ ■ ■ X d d , 
j! = (ji!)(j 2 !)---(jd!), and IP = {dXl)H Xy 2 ... {dXd) u - 

Throughout the paper we assume a 6 lR d , rj G M d and K > 0. We also assume that A 
and B are d x d complex invertible matrices that satisfy 



A 1 B 
A* B 



B* A 
B* A 



0, 
21. 



(2.1) 



These conditions guarantee that both the real and imaginary parts of BA are symmet- 



A A*. 



ric. Furthermore, Re BA~ X is strictly positive definite, and (Re BA 

Our definition of (fj(A, B, h, a, 77, X) is based on the following raising operators that 
are defined for m = 1, 2, . . . , d. 



A m {A,B,h,a,riY 



1 



2h 



d 



Y B nm (X n -a n ) - i 53 A nm (-iU-^zz r) n ) 



71=1 



n=l 



0X n 



The corresponding lowering operators A m (A, B, h, a, rj) are their formal adjoints. 

These operators satisfy commutation relations that lead to the properties of the 
fj(A, B, h, a, rj, X) that we list below. The raising operators A m (A, B, h, a, 77)* for 
m — 1, 2, . . . , d commute with one another, and the lowering operators A m (A, B, h, a, rf) 
commute with one another. However, 

A m (A, B, h, a, rj) A n (A, B, h, a, r/)* — A n (A, B, h, a, rj)* A m (A, B, h, a, rj) = S m>n . 



Definition For the multi-index j = 0, we define the normalized complex Gaussian wave 
packet (modulo the sign of a square root) by 

<p (A, B, h, a, 77, X) = Tr-^/r^detCA))- 1 / 2 

x exp { - ( (X - a), B A' 1 (X - a) ) /(2H) + % (77, (X - a)) /%} . 

Then, for any non-zero multi- index j, we define 

1 



<Pj(A, B, h, a, 77, 



Ai(A, B, h, a, t?)*) Ji ( A 2 (A, B, h, a, 77) 



*\J2 



x (A d (A,B,h,a,7])*) jd ip (A,B,h, a, r/. 
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Properties 1. For A — B — I, % — 1, and a = rj = 0, the <fj(A, B, h, a, rj, • ) are just 
the standard Harmonic oscillator eigenstates with energies \j \ + d/2. 

2. For each admissible A, B, %, a, and r), the set { (pj(A, B, %, a, rj, ■)} is an orthonormal 
basis for L 2 (R d ). 

3. In [|T^], the state <fij{A, B, %, a, rj, X) is defined as a normalization factor times 

Hj(A; h- 1/2 {A]' 1 (X - a)) <p (A, B, h, a, rj, X). 

Here Tlj(A; y) is a recursively defined |j| th order polynomial in y that depends on A only 
through Ua-i where A = \A\ Ua is the polar decomposition of A. 

4. By scaling out the \A\ and h dependence and using Remark 3 above, one can show that 
TCj(A] y) e~ v I 2 is an (unnormalized) eigenstate of the usual Harmonic oscillator with energy 
|i| + d/2. 

5. When the dimension d is 1, the position and momentum uncertainties of the 



(fj(A, B, h, a,r),-) are J (J + 1/2)% \A\ and J(j + 1/2)% respectively. In higher di- 



mensions, they are bounded by y(\j\ + d/2)% \\A\\ and y(\j\ + d/2)% \\B\\, respectively. 

6. When we approximately solve the Schrodinger equation, the choice of the sign of the 
square root in the definition of ipo(A, B, %, a, r), ■) is determined by continuity in t after an 
arbitrary initial choice. 

The following simple but very useful lemma is proven in 

Lemma 2.1 Let P\j\< n denote the projection onto the span of the <Pj(A, B, %, a, T], ■) with 
\j\ < n - 

(X — a) m P\j\< n = P\j\<n+\m\ (X — a) m P\j\<ni (2-2) 

and 

(X-arP mn \\ < (V2Sd|M||) H (<!i±Miy /2 . (2.3) 

In the Born-Oppenheimer approximation, the semiclassical dynamics of the nuclei is 
generated by an effective potential given by a chosen isolated electronic eigenvalue E{X) of 
the electronic hamiltonian h{X), X e lR d . For a given effective potential E(X) we describe 
the semiclassical dynamics of the nuclei by means of the time dependent basis constructed 
as follows: 

By assumption H 1; the potential E : H C lR d — » 1R is smooth and bounded below. 
Associated to E(X), we have the following classical equations of motion: 

a (t) = v(t), 

m = -VE(a(t)), 

A(t) = iB(t), (2.4) 
B(t) = iE®(a(t))A(t), 

S(t) = ^-E(a(t)), 



where E^ denotes the Hessian matrix for E. We always assume the initial conditions A(0), 
B(Q), a(0), 7/(0), and 5(0) = satisfy (|2p . 

The matrices A(t) and -B(t) are related to the linearization of the classical flow through 
the following identities: 

Because E is smooth and bounded below, there exist global solutions to the first two 
equations of the system Q2.4Q for any initial condition if H = lR d . From this, it follows 
immediately that the remaining three equations of the system ( |2.4| ) have global solutions. If 
H 7^ iR d , for any initial conditions, there exists a < T < oo so that solutions to the system 
( j2.4p exist for any time t 6 [0, T], T is finite if and only if the solution a(t) corresponding 
to the chosen initial condition leaves the set S in finite time. 

Furthermore, it is not difficult [ TT| , |T2| to prove that conditions Q2.1| ) are preserved by the 
flow. 

The usefulness of our wave packets stems from the following important property If 
we decompose the potential as 

E(X) = W a (X) + Z a (X) = W a (X) + (E(X)-W a (X)), 

where W a (X) denotes the second order Taylor approximation (with the obvious abuse of 
notation) 

W a {X) = E{a) + E {1 \a) {X - a) + E {2 \a) {X - a) 2 /2 
then for all multi-indices j, 

() r 3 «W/* w(A(t),B(t),h,a(t),r,(t),X) 



ih — 

dt 



= (-y A + W a{t) {X)^j [e^ h ^(t),5(t),M(*),^),*) 

if Ait), Bit), a(t), r](t), and S(t) satisfy ( |2.4| ). In other words, our semiclassical wave packets 
ifj exactly take into account the kinetic energy and quadratic part W a u){X) of the potential 
when propagated by means of the classical flow and its linearization around the classical 
trajectory selected by the initial conditions. 

In the rest of the paper, whenever we write ipj(A(t), B(t), h, a(t), r](t), X), we tacitly 
assume that A(t), B(t), a(t), r](t), and S(t) are solutions to ( |2.4|) with initial conditions sat- 
isfying (PI). 



3 The Born— Oppenheimer Expansion in Powers of e 

In this section we derive an explicit formal expansion in e for the solution to the molecular 
Schrodinger equation by means of a multiple scales analysis. This asymptotic analysis is 
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similar to that performed, e.g., in We discuss this in detail because we need more 

detailed information on the structure of successive terms in the expansion. 

We start with the molecular Schrodinger equation for d nuclear configuration dimensions, 



i e 



~dt 



A x ^ + h{X)V. 



(3.i; 



We consider the isolated, multiplicity one, smooth eigenvalue E(X) of h(X) of hypothesis 
Hi. For the moment we assume E(X) is well defined on all of IR d rather than just on a 
subset H C IR d . Later we introduce a cut-off function to take care of the general case. We 
consider the solution a(t), r)(t), A(t), Bit), and S(t) to the system (|2.4|) of ODE's. Then, 
we choose the phase of the eigenfuction $(X, t) so that 



($(X,t), (i^-+irj(t)V x )^X,t)) 



0. 



(3.2) 



This can always be done. See, e.g., [10 



The multiple scales analysis consists of separating the two length scales that are important 
in the nuclear variable X. The electron wave function is sensitive on an 0(1) scale in this 
variable, so X, or equivalently, X — a(t) is relevant. The quantum mechanical fluctuations 
of nuclear wave function occur on an 0(e) length scale, so (X — a(t))/e is also relevant. 
We replace the variable X by both w = X — a(t) and y = w/e, and consider them as 
independent variables. This leads to the new problem of studying 



1 e 



~dt 



A, 



e 3 V,„ • V, 



A y + ie 2 7](t)-V w + i€T}{t)-V % 



+ [ h(a(t) + w) - E(a(t) +w)) + E(a(t) + ey) 



(3.3) 



We easily check that if ^(w,y,t) solves ( jOD then ^(X - a(t), (X - a(t))/e,t) solves (|3TTD. 
We define $(w,t) = ${X,t). Then (g]|) becomes 



<9 

(Q(w,t), i — $(w,t)} He 



■ 0. 

We seek solutions to ( |3.3|) of the form 

^(w,y,t) = , 
This requires 4>{w,y,t) to satisfy 



(3.4) 



iS(t)A 2 Jv(t)-y/e 



[w,y,t). 



1 e 



dt 



A,. 



^A y + ^E^(a(t))y* 



e 6 V w -V y + 
+ [h(a(t) +w) - E(a{t) +w)\ 

+ ( E(a(t) + ey) - E(a(t)) - e£« (a(t)) ■ y - e 2 E^ (a(t)) 



(3.5) 



2! 
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where here and below we make use of the shorthand notation 

in the usual multi-index notation. We next assume that <p(w,y,t) has an expansion of the 
form 

(p(w,y,t) = (f) (w,y,t) + e0i(io,y,t) + e 2 (j) 2 (w,y,t) + ■■■ 
We further decompose each <p n as 

(p n (w,y,t) = g n (w,y,t)$(w,t) + <f>i(w,y,t), 

by projecting into the $(«;,£) direction and into the orthogonal directions in H. e i- 

We substitute this expansion into ( |3~5D and equate terms of the corresponding powers of 

e. 

Order 0. The zeroth order terms require 

[h(a(t) +w) - E(a(t) +w)\ cj) (w,y,t) = 0. 

This forces 

</>o(u>,y,t) = 0. 
Order 1. The first order terms require 

[h(a(t)+w) - E(a(t)+w)) (p!(w,y,t) = 0. 

This forces 

<f>i{w,y,t) = 0. 
Order 2. The second order terms require 



= ( --A, ; ^E^(a(1))^ ) o n +- [/,(„(,) + ir ) £(,(/) „-)],;,. 

We separately examine the components of this equation in the $ direction and in the or- 
thogonal directions. By (|3.4j), this yields the two conditions 



and 



We arbitrarily choose go to be the following ^-independent particular solution of ( |3.6| ): 

g (w,y,t) = e- d ' 2 coj <Pj(A(t), B(t),l,0,0,y), (3.8) 

lil<J 
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^ = (-iA, + ^)(a(t))|), , (3.6) 



<9$ 

[ h(a(t) +w) - E{a{t) +w)}<f) 2 = igo^- (3.7) 



where cqj = Cj is determined by the initial conditions. 

We let the Hilbert space 7i^j be the subspace of TC e \ orthogonal to t). The restriction 
of [h(a(t) + w) — E(a(t) + w) ] to Ti^x is invertible, and we denote the inverse by r(w, t) = 
[h(a(t) + w) — E(a(t) + w)}' 1 . With this notation, equation ( |3.7|) forces 

02 ( w > y>t) = 1 9o(w, y, t) r(w, t) -^{w, t) 

= e- d / 2 ]T d 2J (w,t) <pj(A(t), B(t),l,0,0,y), (3-9) 

\j\<J 

where 

d 2 ,j(w,t) = c 0J r(w,t) — (w,t) (3.10) 
is Hgi^valued. 

Order 3. The third order terms require 

3 

-V w - V v 0o + ^ 3 )(a(t))|-0o + [h(a(t)+w) - E{a{t) + w)} 3 - 

We separately examine the components of this equation in the $ direction and in the or- 
thogonal directions. By (|3.4|), this yields the two conditions 



y 3 



(V^o)-(^ V,$) + ^ 3 )(a(t))|- po, (3-11) 



and 



9$ 

[%(*)+«;) - £(a(t)+w)] 3 = ^i^- + (V^ ) ■ (PlV^), (3.12) 

where P±(w,t) is the projection in H e z onto 

The solution to ( BTTTj ) with gi(w, y, 0) = can be written as 

9i(w,y,t) = e~ d ' 2 ]T c ld {w,t) <Pj(A(t), B{t),l, 0,0, y), 

\j\<J+3 

for some coefficients Cij(w,t). 
Equation ( |3.12|) determines 

03 = r(w,i) — (w,t) + {V y go)(w,y,t) ■ (P±(w,t)V w $(w,t)) I 



e 



-rf/2 



£ d 3 j(w,t) V3 (A(t),B(t),l,0,0,y), 



|il<-/+3 
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where 

d 3 ,j(w,t) = i (r(w,t)$(w,t)j c ltj (w,t) 



\9\<J 



d 

Here and below = — . 

dt 



Order n. The n order terms require 



-A y + £ (2) (a(t))^- n _ 2 - -A^„_4 - V w -V y (f) n ^ 



dt \ 2 u v v " 2! J r ™ 2 

n m 
+ J] # M («(*)) ^7 0n-m + - E(a(t)+W)] (j>, 



m=3 

The components of this equation in the Q(w,t) direction require 

= - - A w g n - 4 - ( $, V w $ ) • (V w g n - 4 ) - - ( ) #„_ 4 



m=3 m - 

-I($,A^_ 4 > - ($, V u -V^_ 3 > + i(^,0i_ a >. (3.13) 

8(f) 1 - <9$ 

Note that the last term has been transformed from — i ( $, — ) to i ( — — , dab 9 ). The 

dt x dt n ~ 2 1 

equivalence of these expressions follows from differentiation of ( $, <^_ 2 ) = with respect 
to t. 

The components orthogonal to $(iu,t) require 
[ /i(a(t) + w) - E(a(t) + w) } (j) n 



m=3 

+ -P ± A TO 0^_ 4 + (Pj_ V w $) ■ (V u fl n -4) + \ (Pj. A w $) <? n -4 

0n-3 + (fiV^)-^) + *-fc9n-2- (3.14) 

Equation ( |3.14j ) determines 4>^(w, y, t) by an application of [h(a(t) + w) — E(a(t) + w)]~\ 
It is easily checked that the solution to ( |3.13|) with g n _ 2 (w,y,0) = has the form 

g n -2{w,y,t) = e- d l 2 E c n _ 2j >,t) ^(A(t),B(t),l,0,0,y), (3.15) 

|i|<J+3n-6 
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for some coefficients c n _2j(w, t), and that the y dependence of the vector 0^ has the same 
form, with other coefficients depending on (w,t), i.e., 



<f£(w,y,t) = e- d/2 Y d n>j {w,t) ^(A(t), B(t), 1,0,0, y), 



(3.16) 



\j\<J+3n-6 



where the d n j(w,t) take their values in the electronic Hilbert space. 

Equations (|3.13|) and fl3.14| ) determine c n _2,j and d n j. When recursively solving these 
equations, we must determine d n j before c n j because the right hand side of ( |3.13| ) (with 
n — 2 replaced by n) contains 0^. 

The solution to ( |3.14|) in terms of the d n j, is 



d nJ (w,t) = Y Ai(w,t), 



(3.17) 



where 



Ai(ty,t 
A 2 {w,t 

A 3 (w,t 
A±(w, t 
A 5 (w,t 
A e (w,t 

A 7 (w,t 

A 8 {w,t 



ir(w,t) P±(w,t) d n - 2tj (w,t) 
(D m E)(a(t)) 



- E 

3<\m\<n 



ml 



Y V m Vq) r(w,t) d n -\ m i q (w,t) 

q\<J+3(n-\m\-2) 



-r(w,t)P±(w,t) (A w d n ^ tj )(w,t) 
r(w,t) P±(w,t) (V w $) ■ (V w c n - 4:j )(w,t) 



-r(w,t)P X (w,t) (A w $) Cn_4j(tW,*) 

Y r(w,t) Pj_(w,t) (<fj, Vyifq) {V w d n - 3t q)(w,t) 

\q\<J+3(n-B) 

Y r(w,t) P±(lV,t) (V w $) (<pj, S/yVq) Cn-3, q (w,t) 
\q\<J+3(n-3) 

ir(w,t)Pj_(w,t) $(w,t) c n _ 2J (w,t). 



Similarly, the solution to ( |3.13j ) in terms of the c nj - is obtained by integration with respect 
to t of i c n j(w, t), where 



ic nJ (w,t) = Y r ;(?M) 

8=1 



where 



r 2 (w,t) 
r 3 ( w ,t) 



(3.18) 



1 



(A w c n - 2 ,j)(w,t) 



($, V w $) • (V w c n _ 2J ){w,t) 

\ ($, A w $) Cn- 2 j(w,t) 



12 



|g|<J+3(ra-l) 

T 5 (w,t) = - ]T ( $ ' V » $ ) ' (<Pj> V ^<?) Cn-l,q(w,t) 

|g|<J+3(n-l) 

i v ^ ^ (L> m £)(a(t)) _ , . , 

3<M<n+2 | 9 |<J+3(n+2-m) m ' 

r 7 (w,t) = l - ($, (A M d n _ 2ii )(w,t)) 

\q\<J+3(n-3) 

T 9 (w,t) = i ($, d nJ (w,t)). 

4 The Main Result 

We introduce a C°° real valued cut-off function F : lR d — ► iR that equals 1 in a neighborhood 
of the origin and equals zero away from the origin. More precisely, we choose < b Q < b\ < 
oo, such that 

supp (d w .F)(w) C {w G lR d : b < \w\ < &i}, 

for any i G {1, . . . , d}, and such that for any tE.fl, all quantities appearing in the above 
expansion are well defined for w G lR d with \w\ < b\. Here f2 is a particular simply connected 
open complex neighborhood of the real interval [0, T] that we construct in Section [5] under 
hypotheses H and Hi. 



We define our approximate solution to (|3.1|) at order iV by the following expression: 



* N (w,y,t) (4.1) 

/ N N+2 \ 

= F(w)e* s ^ 2 e^t>y/t ^ e n g n {w,y,t) <$>{w,t) + £ e n ^(w,y,t)\. 

\n=0 n=2 / 

We prove in Section |7.2| that this quantity agrees with an exact solution up to an error whose 
norm is bounded by for t G [0, T]. 

We emphasize that once the molecular hamiltonian h(X) and its spectral data E(X), 
$(X) are given, the only arbitrary input of the above derived expansion consists of the set 
of coefficients c j, \j\ < J. We note that at time t — 0, we have c n j(0,w) = for all n > 1. 
Thus, at t — 0, the approximation reduces to 

V N (w,V, 0) 

/ iV+2 \ 

= F{w) e iS ^ e^°>^ g o (0,y,0)$(w,0) + £ e n ^(w,y,0) . 

\ n=2 / 

This expression is completely determined by go{0, y, 0), the nuclear part of the wave function 
parallel to the chosen electronic level at time 0. 
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As is usual in the study of adiabatic problems, in order to get accurate information on 
the evolution of an initial wave function associated with a specific electronic level, one needs 
to include a higher order component perpendicular to that electronic level. This higher order 
part is completely determined by the parallel part. Here it is given (up to phase and cut-off 

iV+2 

functions) at time by e n <j)^(w,y,0). 

n=2 

We now state our main Theorem: 



Theorem 4.1 Assume hypotheses Ho and Hi and consider the above construction. For 
all sufficiently small choices of g > 0, there exist C{g) > and T(g) > such that, for 
N(e) = [g 2 /e 2 ], the vector ty*(X,t,e) = ^ N(e) (X - a(t), (X - a(t))/e, t) satisfies 



L 2 (R d ,n el ) 



< C{g) e 



-r( 9 )A 2 



for allt G [0, T], as e -»• 0. 

Moreover, we have the following exponential localization result. For any b > and a suf- 
ficiently small choice of g > (that depends on b), there exist c(g) and 7(g) > 0, such 
that 

\ 1/2 



/ 

J \x—a( 



i\x-a(t)\>b 

for allt G [0, T], as e -> 0. 



\V*(X,t,e)\\ 2 Hei dx 



< c{g) e" 7(9)/e , 



The strategy of the proof is as follows: We consider the approximation ^! ^{X — a(t), (X — 
a(t))/e, t) and the exact solution to the Schrodinger equation with the same initial conditions. 
We estimate the norm of the error (that is the difference between these two quantities) as 
a function of both iV and e. Apart from some subtleties, the norm of the error is bounded 
by Ce N \tN 1 / 2 ) n , for some constants C and r > 0. We minimize the error estimate over all 
choices of N. This yields iV ~ g 2 /e 2 , for sufficiently small g > 0, and an estimate of order 
e -r(s)A f or tfig norm of the error. 

We prove two extensions of this result in Section || In the first extension, we consider the 
validity of our approximation on the Ehrenfest time scale, i.e., when T = T(e) ~ ln(l/e). In 
the second extension, we study the dependence of our construction on J, in order to extend 
our main result to a wider class of initial conditions. We refer the reader to Section |] for 
the precise statements. 



5 Analyticity Properties 

Our estimates depend on analyticity in t G £1 of the vectors c n (w,t) G Z 2 (N d ,(Z7) and 
d n (w,t) G l 2 (N d ,H e i), where f2 is the particular simply connected open complex neigh- 
borhood of the real interval [0, T] mentioned at the beginning of Section |j. 

To construct Q, we begin with several observations. Our hypotheses imply that the 
eigenvalue E(X) is analytic in Ey, so the solutions a(t), r](t), A(t), B(t), and S(t) are well 
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defined for all t G [0, T]. Moreover, by standard arguments ||, these functions all have 
analytic continuations from [0, T] to a simply connected open set that contains [0, T]. 
We assume without loss of generality that Oi = Jli, where fli denotes the conjugate of Qi. 

We note that A*(t) and B*(t) also have analytic continuations from [0, T] to flj. To see 
this for A*(t), note that for t G [0, T], A*(t) = A*(t), and A*(t) has an analytic continutation 
to Q\. The argument for B*(t) is similar. 

It now follows easily from the definitions that for each X, ipj(A(t), B(t), e 2 , a(t),rj(t),X) 
and ipj(A(t), B(t), e 2 , a(t), rj(t), X) have analytic continuations from [0, T] to some simply 
connected open set f2 2 . For t G [0, T], the real part of is strictly positive. This 

positivity will remain true for the real part of the analytic continuation of B(t)A(t)~~ l on 
some simply connected subset Q C Oi fl ^2 that contains [0, T]. We assume without loss 
of generality that Q = Q and we can assume that Q has the form {t : —a < Ret < 
b and |Imt| < c} where a > c > and b > T + c. It follows that for £ £ f2, both 
(fj(A(t), B(t), e 2 , a(t), rj(t), x) and (fj(A(t), B(t), e 2 , a(t),r](t), x) have analytic continuations 
from [0, T] to as elements of L 2 (lR d ). 

Using these results and carefully examining the constructions of the vectors c n (w,t) and 
d n (w,t), we see that they are analytic in t for t G Q, also. 

Our hypotheses on h(-) and the above results also show that each of the following quan- 
tities is analytic in t for t G f2 and each fixed w G S«5 C for sufficiently small 5: 

r(w,t) = [h{a{t)+w) - E(a(t) + w)]' 1 , 

$(w,t), 

(D^){w,t), for|a|<2, 
D%E(a{t)), for all a 
P±(w,t). 

By explicit computation of the phase corresponding to (|3.2|) it is easy to check that <&(w,t) 
and its derivatives are also analytic for t G O. 

Moreover, if fi(w, t), (i in some finite set) represents any of these quantities, fa is analytic 
in w G Ej, for any fixed t G £1. Thus, by the Cauchy integral formula, we can assume that 
the following bounds hold (with the appropriate norm in each case): 

ii wok oil < ^r 1 

for some q, Gj, w G Tig, an d a ranges over all multi- indices. 
We can assume here that all Gi < D2 for some constant 

D 2 > 1, 

and we associate the prefactors Cj in (|5.1|) with the different functions according to the rules 

Ci <r-¥ rP± C2 <r-> $ 

c 3 <-> $ c 4 <-> V w $ 

c 5 <-> A w $ c 6 <-> 

4 <-> V M $) c' 5 ^ ($, A^$). 
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6 Structure and Estimates of the c n (w,t) and d n (w,t) 

In this section, we decompose the functions g n and <fi^ of Section § into pieces, each of which 
satisfies various estimates. 

Throughout this section, all ^-dependent quantities are defined for w in the support of 
the cut-off function F. Furthermore, all the results of this section are claimed to hold only 
on the support of F. 

Our decompositions of g n (w,y,t) and (f>^(w,y,t) have the following forms: 

9n(w,y,t) (6.1) 
= £~ d/2 E E E E c n , PiiW >,i) <Pj(A(t), B(t),l,0,0,y). 

$eB„,i p<n \i\+k<p+f \j\<j +n +2(p-\l\-k) 

and 

<tt(w,y,t) = (6.2) 
e - d /2 E E d niPMj (w,t) <Pj(A(t), B(t),l, 0,0 )2/ ). 

/3eB n , 2 p<n-l \i\ +k < p+ S-z± \j\<J+( n -l)+2(p-\l\-k) 



In (5TT), n, and p are non- negative integers; j and I are multi- indices; and the index (3 
runs over a finite set B n> \. The number J is fixed by the initial conditions. Each c^pj^pj is 
a complex valued function. 

In (|6.2j ), n > 2, k and p are non-negative integers; j and / are multi- indices; and the index 
(3 runs over a finite set £y 2 - Each dn,p,i,k,pj{ w ,t) takes values in 7i e /. 

We let c n:P: i : k,/3(u!,t) and d n ^k,p{w,t) respectively denote vectors in l 2 (N d ,(U) and 
l 2 (N d ,H e i) whose components are c ntPt i : k t pj(w,t) and d ntP j : k,p,j(w,t). 



The crucial step in the proof of Theorem |4 . X| is the following: 

Proposition 6.1 There is a recursive construction of the coefficients Cn, p ,i,k,p,j(w->~k) and 

d n , P ,i,k,/3,j(w,t) form on the support of F. 

The indices for c n ^i^,p,j{w ,t) are non-negative and satisfy 

(3 e B n>1 
p < n, 

\l\ + k < p+-, 

\j\ < J + n + 2(p—\l\ — k). 
The indices for d ntP j t k t /3j(w,t) are non-negative and satisfy 



n 


> 


2 


P 


e 




P 


< 


n — 1 


l\ + k 


< 


n 

p H 


\J\ 


< 


J+(? 



2(p- |/| - k). 
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Moreover, the following conditions are satisfied: 

i) For any n > 0, c nfl> i jk)/3>j (w,t) = 0. 

ii) There exists Kq > 0, such that the number of terms in both of the sums (6A) and (6.i) 
is bounded by e K ° n . 

Hi) For t G Q, let dist(i) be the distance from t to the complement of Q. The coefficients 
c n,p,i,k,p[ w 1 1) and dn,p,i,k,p( w )t) are analytic for t e Q, and there exist constants D\ and D 2 , 
such that 



\\{Dlc n ^i A(5 ){w,t)\\ 

a + l)\ \t\ p k k 



(6.3) 



< D 1 D 



H + |Z|+4n 



[l + \a\) d + l p\ dist(t)' 



{J + n + 2(p- \l\ - k))\ 
J! 



1/2 



and 



\\{Dy n>Pil>k)l3 )(w, 



(6.4) 



|a|+|J|+4(n-l) (a + l)\ \t\ p 



[l + \a\) d+1 p\ dist(t)' 



(J+(ra-l) + 2(p-|Z|-fc))! 
J! 



1/2 



Remark The complicated estimates ( |6.3|) and (6.4) are motivated by estimates used in 
semiclassical approximations and adiabatic approximations. The factors on the right hand 
sides that explicitly involve J, n, and p occur in the semiclassical paper [|I~5 |. The factors 
that involve a and I appear in the adiabatic paper f25|. The factors that involve k occur in a 



proof of the adiabatic results of |55| that are based on Cauchy estimates instead of Nenciu's 
lemma p5f (that we generalize below as Lemma |6.4j ). We were unable to prove Proposition 
|6.1| without using a combination of all of these techniques. We estimate adiabatic error terms 
by using Nenciu's approach in the w variable and Cauchy estimates in the t variable. 



6.1 The Toolbox 



To prove Proposition |6.1| , we repeatedly use the following very handy lemmas, whose proofs 
are given in Section [| The first two lemmas deal with basic properties of analytic functions 
of one variable and are consequences of the Cauchy integral formula. 



Lemma 6.1 For k = 0, define k k = 1. Suppose g is an analytic vector-valued function on 
the strip Ss = {t : |Imt| < 5}. If g satisfies 

\\g(t)\\ < Ck k (5 - |Imt|)- fc , 

for some k > 0, then g' satisfies 

\\g'{t)\\ < C(k + l) k+1 (6 - \lmt\)- k -\ 

for all t e Ss- 
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Lemma |6J] has a generalization to regions other than infinite strips. The generalization 
is needed if one wishes to study problems where analyticity holds only in a neighborhood of 



a finite time interval. The proof of the generalized lemma is similar to that of Lemma O 
but involves slightly more complicated geometry. The precise statement is the following: 



Lemma 6.2 For k = 0, define k k = 1. Suppose g is an analytic vector-valued function in 
an open region Q C (D . For t G Q, let dist(t) be the distance from t to Q c , the complement 
of ft. If g satisfies 

\\g(t)\\ < Ck k (dist(t))- k , 

for allt G f2 and some k > 0, then g' satisfies 

\W(t)\\ < C(k + l) h+1 (dist(t))- k -\ 

for all t G Q. 

The next lemma gives estimates on indefinite integrals of certain analytic functions under 
stronger assumptions on the domain Q. 



Lemma 6.3 Suppose f is an analytic vector-valued function in an open region C (U. For 
t G Q, let dist(t) be the distance from t to Q c . We assume the domain is star-shaped with 
respect to the origin and that the origin is the most distant point to fl c , i.e., dist(O) > dist(t), 
for all t G Q. Moreover, we assume that dist(t) is monotone decreasing along any line 
emanating from the origin. If f satisfies 

\\f(t)\\ < C \t\v (dist(t))- fe , 



for all t E Vl and some k > 0, then J * f(s) ds satisfies 



f{s)ds 



< c LL_ (dist(t))- fc , 



for all t G O. 



Remark: In our situation, examples of sets f2 we can use that satisfy the conditions of 



Lemma |6.3| are infinite symmetrical horizontal strips or the rectangular regions chosen in 
Section [5]. 



A fourth tool we repeatedly use below is a multidimensional generalization of a lemma 
used in ||25|| . We warn the reader that the symbol for a norm means different things in 
different contexts, e.g., for scalar-valued, operator-valued, and vector-valued functions, it 
respectively means absolute value, operator norm, and vector space norm. 
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Lemma 6.4 The quantity 

" = S " P (1 + Wr aT^r (TT^ (6 - 5) 

is finite. 

Let £ 6e an open subset of(U d . Suppose M(-) G is scalar-valued or operator-valued, 

andN(-) G C°°(S) zs either operator-valued or vector-valued. Assume these functions satisfy 

\\ (D<* M) (x) \\ < m(x)a(x)^ ( i+j^L (6.6) 

HD°N)(z)\\ < n{x)a{xt+«\ (/"^L (6-7) 
/or i£E, a// multi-indices a, and some fixed multi-indices p and q. Then 

|| (D a (MN))(x) || < m(x) n(x) v a{x)\ a+p ^ ^4^M (6-8) 
for each multi-index a, where v is defined by (W 



6.2 Proof of Proposition 

We prove Proposition by induction and begin with the case n = 0. We construct 
co,o,o,o,/9j = c o,j with (3 = 1 G £>o,i = {1}. We note that there is no d ntP) i t k,/3(w,t) for n < 1; 
the inequalities for its indices in the conclusion to the proposition cannot be satisfied by 
non-negative integers. Whenever d ntPt i t k,/3{w,t) with n < 1 appears in any of the formal 
calculations below, it is understood to be zero. 

We now assume that the estimates ( |6.3|) and (|6.4|) on Cm, p ,i,k,p{w ,t) and dm^j^^iw.t) are 
true for all m < n — 1 and prove they still hold for m — n. 

Our strategy is to show that each contribution Aj and Tj consists of a finite sum of terms 
that satisfy the required estimate. We estimate the number of terms by a separate argument. 
Our main tools are Lemmas 12. 1L 16. 2L B731 and 



The index (3 must be considered when counting the number of terms, but it plays no 
role in the estimates of the individual terms. To simplify the notation, we drop it while 
estimating the terms. 

The Term A x 



We begin by considering the contribution to (|6.4j ) from the term Ai in (|3.17 ). 

By induction, each d n - 2 , P ,i,k,i3( w -> t) is analytic for t G O and has a p th order zero at t — 0. 
It follows that d n -2, p ,i,k,p(w,t) = t p fit), where / is analytic in Q. When we take the time 
derivative, we obtain two terms, pt^ 1 fit) and t p f(t). These, respectively, give rise to two 
terms d n!P -i,i, k ,p>{w,t) and d n ^ k+1 ^>{w,t). 
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We consider all ^-derivatives of Ai(w,t). We apply the induction hypothesis, Lemma 



672] , and Lemma |6.4| to obtain 

DZ (rP ± (w,t)d n _ 2 ^i,k(w,t)^ 

|a|+|i|+4(n-3) + Q! 



'(J+(n-3) + 2(p-\l\-k))\ 



;i + H) d+1 



J! 



t p (Jfe + l 



,fc+i 



(p- 1)! dist(t) fc p! dist(t) 



fc+i 



(a + Z)! 



|t| p 



1 + p'! dist(t)' 



+ c 1 vD 1 D~ 8 D 



\a\+\l\+4(n-t) (a + Z)! |t| 



k' 



I {J + {n - 1) + 2{p' - \l\ - k))\ 

; (J+ (rc- 1) + 2(p- |Z| - fcQjT 
J! 



(1 + p! dist(t) fc ' V 

with p' = p — 1 and fc' = k + 1. We check that 

p < n — 3 < 7i—l, 

p' < n — 4 < n — 1, 

|Z| + fc < p+(ti-3)/2 = p' + (n-l)/2, 

+ < p+(ti-3)/2 + 1 = p+(n-l)/2, 

and the ranges of the components of each vector satisfy 

\j\ < j + {n-l)+2{p' -\l\-k), 

\j\ < J + (n-l)+2(p-\l\-k'), 
as required. Hence, we get the desired bound for each of the two contributions, provided 



Dl > a v. 



The Term A 2 



In the analysis of this term, we encounter an infinite matrix that represents multiplication 
by y m in the basis of semiclassical wave packets. We denote this matrix by (<p, y m ^p). Its 
entries are (ipj, y m (f q )(t), for multi-indices m,j,q G N d . We recall that Lemma |2~T1 gives 
bounds for these matrix elements and also states that (<pj, y m (p q )(t) = if | \j\ — \q\ | > \m\. 

We adopt the analogous notation for the infinite matrix (<p, D™(p) that represents the 
operator D™ in the basis of semiclassical wave packets. 

We define do = V2d. Then, using ( |5.1|) , Lemmas pTTj |673| , |B~2| , and|0], and some algebra, 
we obtain 

E <V- irs)U-Pi.)UrJ),h.- l , l ,.n,Uc.t) 



m=3 \m\=m 



ml 
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< E (1 + f „)« m , 

m=ci m|=m v ' 



Ceciz/D^m! /( J + (n - 1) + 2(p - |Z| - fc))! 



J! 



|a]+|i]+4(n-l-m) (<* + /)! |t| P & A 



(1 + \a\) d+1 p\ dist(t) fc 



< 2^ ^lC 6 CiZ/ — i>» 2 



m=3 |m|: 



(1 + \a\) d+l p\ dist(t) fc 
/ (J + n -l + 2(p- |l| - fc))[ 

We also verify the constraints on the parameters and components of the vectors: 

p < n — 1 — rh < n — 1, 
|Z| + fc < p + (n-l-m)/2 < p + (n-l)/2, 

|j| < J+(n-m-l) + 2(p-|Z|-ife)+m < J + (n - 1) + 2(p - \l\ - k). 

Hence, we see that each contribution from A 2 (w,t) satisfies the required bound, provided 
the following two conditions are fulfilled 

Dl > do\\A\\, 

Dl > {d Q \\A\\fc, Cl u. 

There are 



E E 1 < E 1 

m=3 \m\=fh \m\<n 



n + d 

d 



< a e° 



such contributions, where a > can be chosen arbitrarily small (see [H 
The Term A 3 

For this term we make the laplacian explicit and write 



A„A_ 4jMfc O,t) = ( D lidn-4,p,l,k)(w,t). 



i=l 



We introduce = I + (0, 0, . . . 0, 2, 0, . . . , 0), where the 2 sits in the ith column. 
We then estimate 



DZ - (rP ± (w,t)A w d n - 4 ,p,i,k( w > t )) 



^ 1 ^ n n H + |Z|+2+4(n-5) (a + h,2V- \t\ P k k 

^ O E C l V D ± D 2 



1=1 



{l + \a\) d+1 p\ dist(t) ? 



(J + (n — 5) + 2(p - \l\ - k))\ 



J\ 
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* C 1 UD 1 |a|+|i il2 |+4(n-l) (<* + ^2)! W 

- 2^ o nie ■ 



, 2L> 2 16 2 (l + H)^ 1 p! dist(t) fc 



x 



(J + (n — l) + 2(p - |Z i)2 | - fc))! 



J! 



Again, the constraints are satisfied since 

p < n — 5 < n — 1 

< p+(n-5)/2 + 2 = p+(ra-l)/2 
|j| < «/ + (n — 5) + 2(p — |Z| — k) = J + (n — l) + 2(p— \l ij2 \ - k) 

and each of the d contributions stemming from A 3 (w,t) satisfies the required estimate, 
provided 

> Cl u/2. 

We estimate each of the remaining terms Aj(w, t), i — 4, . . . , 8, in the same fashion, using 
the same tools. Since this is straightforward, we only outline the arguments. 

The Term A 4 

We expand the dot product 

d 

i=i 

and use the definition 

l iA = l+ (0,0,... 0,1,0,..., 0), 

where the 1 sits at the ith column. Recall that the estimates on the c^p^'s differ from 
those on the d mtPt i t kS by a shift of 1 in the m dependence. We have 

\\Dl (rP ± (w,t) V w $(w,t)- ( w ,t)) II 

jLv c 1 c A u 2 D 1 H+|; M | +4 (n-i) + |t| p /c fc 



i=i 



Df 1 (l + lal)^ 1 p\ dist(t) fc 



X 



'(J+(n-l) + 2(p-\l i!l \-k))\ 



with all constraints on \j\,p, |/i,i|, A; satisfied. Thus each of the d contributions stemming 
from A 4 (u>,£) satisfies the required estimate, provided 

Df > cic 4 z/ 2 . 



The Term A 5 
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This term is similar to the previous one. We obtain 

K ( \ rP±(v>, t) (w, t) c n ^ pM (w, t)j 



< 



dC 5 U 2 Di | a |+|i]+4(n-l) (<* + /)! \t\ P k h 



'2m 2 ~° 2 



[l + \a\) d+1 p\ dist(t)' 



x 



(J + (n — 1) + 2(p - \l\ - fc))! 
J! 



with all constraints on \j\,p, \l\,k satisfied. Thus the contribution stemming from A5(w,t) 
satisfies the required estimate, provided 

D 12 > Cl c 5 z/ 2 /2. 



The Term A 6 

At this point the matrices {ip, D yi tp) play a role that we control by the momentum space 
analog of Lemma pTT] . Expanding the dot product and introducing the matrices (ip, D yi ip), 
% = 1, . . . , d we have the following estimate for this term: 



D Z J2 (rP±)(w,t) (ip, D y .ip) D Wi d n - 3jP; i >k (w,t) 



i=i 



< d Q c 1 u\\B\\D 1 ^| a |+|; i|1 |+4(n-i) + \t\ p k k 



i=l 



D 12 



[1 + \a\) d+1 p\ dist(t) fc 



\J + (n-l) + 2(p - - k))\ 



V J! 

with all constraints on |j|,p, satisfied. Thus, each of the d contributions stemming 

from A e (w,t) satisfies the required estimate, provided 

D\ 2 > docii/pll- 



The Term A 7 



Similarly, 



K E (rP ± )(w,t) (cp, D y .<p) (D w .$)(w,t) c n ^ hk {w,t) 



i=l 



< ^0 Cq C 4 ^ 2 ||-B|| £>i ^|a| + |Z|+4(n-l) (« + 

i=l 



1 + p! dist(t)* 



x 



(J + (n — 1) + 2(p - \l\ - k))\ 
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with all constraints on \j\,p, \l\, k satisfied. Thus, each of the d contributions stemming from 
A 7 (w,t) satisfies the required estimate, provided 

D 12 > d dc 4 u 2 \\B\\. 



The Term A f 



Finally, 



K (rP L {w, t) t) c„_ 2 ^ ifc (w, t) 



< ci c 3 v 2 Pi | Q |+|z|+4(«-i) (a + Q! \t\ p k k 

Dl 2 (l + H)^ 1 p\ dist(t) fc 



'(J + (n-1) +2(p-|/|- fc))! 



J! 



with all constraints on satisfied. Thus the contribution stemming from A§(w,t) 

satisfies the required estimate, provided 

D 2 > c\ C3 v 2 . 

We now perform a similar analysis for the quantities Ti(w, t) that appear in the expression 
for c n ^k(w, t). We integrate these terms with respect to t and apply Lemma According 
to the lemma, integration of a term with a given value of p gives rise to a term with p' — p+1 
in the estimates. We also note that the estimates we want to prove for the c's differ from 
those for the d's by the replacement of n — 1 by n. 

The Term / * I\ 

We use the same techniques above to obtain 
ft 1 

D® / - (A w c n _ 2 ,p,(,i)(M), s) ds 
Jo I 



< E 



Dl ^|a| + |i,, 2 |+4n (a + l i)2 )\ t P ' 



. , 2D* ~ A (1 + p'\ dist(t) ? 

We check that the constraints are satisfied: 



' {J + n + 2(jpf-\l it2 \-k))\ 



p' < n — 1 < n, 
\k,2\ + k < p+(n-2)/2 + 2 = p' + n/2 

\j\ < J + {n-2)+2(p-\l\-k) = J + n + 2(p' — \l i>2 \ — k). 

Thus, each of the d contributions stemming from Ti(w,t) satisfies the required estimate 
provided 

D 8 2 > 1/2. 



The Term / * T 2 
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Similarly, with p' — p + 1, 



J ° i=l 



< ^ c^Di £> |a| + | f<[1 | +4 n (a + W W P ' kk J{J + n + 2(p>- - k))\ 



(1 + H)^ 1 p'! dist^)^ 



J! 



with all constraints on satisfied. Thus each of the d contributions stemming 

from T 2 (w,t) satisfies the required estimate, provided 



The Term /„* T 3 

Again, with p' = p + 1 
n 1 



DZ I -($, (A w ®))(w,s) c n - 2 .p.iA w i s ) ds 
o 2 



< c 5 uDl D H+\i\+*n (a + 1)1 \t\ p ' k k J (J + n + 2(p' - \l\ - k))\ 



2DI 



{l + \a\) d+1 p'\ dist(t)' 



with all constraints on \l\,k satisfied. Thus, the contribution stemming from T 3 (w,t) 

satisfies the required estimate, provided 

Dl > c' 5 u/2. 



The Term /„* T 4 

Recall that the matrices (ip, D Vi ip) are controlled by an analog of Lemma |2]T 

rt d 

D Z / Yl D ViV) D Wi c n - ltP ^ k {w,s)ds 

Jo i=l 



1 d \\B\\ Pi nH+ K ll |+4„ (a + k,x)\ \t\* k k j(J + n + 2{p l — — k))\ 



{l + \a\) d+1 p'\ dist(t)' 



J! 



with all constraints on |/i,i|,A; satisfied. Thus each of the d contributions stemming 

from T4(w,t) satisfies the required estimate, provided 



Dl > do\\B\\. 



The Term / * T 5 
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For this term we obtain 

rt d 



K I E ( $ > (D Wi $)(w, s)) (if, Dy.(p) c„_i, Mfc (u>, s) ds 
Jo i=i 



< ^ c^jp ||ff|| gi D \ a \ + \i\ + 4n (a + iy. W_ k* 

1=1 



'(J + n + 2(p' - \l\ - fc))! 



[l + \a\) d+l p'\ dist(t) fc 



J! 



with all constraints on |/|, k satisfied. Thus each of the d contributions stemming from 

F 5 (w,t) satisfies the required estimate, provided 



Do > c' 4 isd \\B\ 



The Term ft T 6 



In this term, we encounter the sum over all previous c's. As in the similar contribution 
from A 2 , we obtain 

K E E f ^§r^ v m v>) ^-^(w, 8 ) 

~ q I I - J lib . 



fh=3 \m\=m 



< y: e d i c * d 

rh=3 \m\=m 



2 pirn 



D 



\a\ + \l\+4n 



{a + l)\ \t\P' k h 



[1 + H)^ 1 p'\ dist(t) fc 



x 



[J + n + 2{p' - \l\ - k))\ 



J! 



We check that the constraints on the parameters and components of the vectors are satisfied 

p' < n — m + 3 < n 
\l\ + k < p+{n-m + 2)/2 < p' + n/2 

\j\ < J + n + 2 + 2(p - \l\ - k) = J + n + 2(p' - \l\ - k). 

Hence we see that each contribution from A 2 (w, t) satisfies the required bound, provided the 
following two conditions are fulfilled 

D\ > d \\A\\ 

D 2 > (d \\A\\fc 6 . 

n 

There are ^ ^ 1 < °"o e CTn such contributions, where a > can be chosen arbitrarily 

rh=3 \m\=m 

small. 



The Term / * T 7 
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This terms depends on the ofs. Recall the estimates are a little different for them. 

d 



K I* Y,\ (H^,s), (D 2 w d n _ 2 ^ k )(w,s))ds 
Jo i=1 2 



< C 2 vD X D \a\ + \l li2 \ + 4n (« + lift)*- W k k 

i=l 



2D 12 2 



;i + \a\) d+1 p'\ dist(t)* 



'(J + n + 2{pf-\l it2 \-k))\ 
J! 



with all constraints on l^l,^ satisfied. Thus each of the d contributions stemming 

from T 7 (w,t) satisfies the required estimate, provided 



D\ > c 2 u/2. 



The Term / * T 8 
Similarly, 

rt d 

J ° i=l 



< y, C 2 Z/rf ||g|| n | a | + |« l , 1 |+4n (a + Z,,i)! \t\ p ' k k 

i=l 



[1 + \a\) d+1 p'\ dist(t)' 



(J + n + 2(j/- \l itl \ -k))\ 
J\ 



with all constraints on \j\,p', \h,x\,k satisfied. Thus, each of the d contributions stemming 
from T 8 (w,t) satisfies the required estimate provided 

D s 2 > c 2 ud \\B\\. 



The Term / * T 9 
Finally, 

D® I (®(w,s), d n7 p 7 i, k (w,s)) ds 
Jo 

< c ^ Dl nM+UI+4« (« + ! ]*K k k l (j + n + 2{pf-\l\-k))\ 
Dl 2 (1 + H)^ 1 p'\ dist(t) fc V J! 

with all constraints on |/|, k satisfied. Thus, the contribution stemming from T 9 (w,t) 

satisfies the required estimate, provided 

D\ > c z v. 

By choosing D 2 large enough, all conditions are satisfied. This completes the induction 
for part iii) of Proposition |6.1| . 

The integration required to construct the c's shows that we obtain non-zero results for 
c n,p,i,k,p f° r n > only when p > 1. This proves part i) of Proposition |6.1| . 

We now turn to the proof of part ii) of Proposition |Q| . 
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6.3 Counting the Number of Terms that Occur in Our Expansion 

In our Born-Oppenheimer expansion, the n th order term has the form 
4> n (w,y,t) = g n (w,y,t)$(w,t)+4>n(w,y,t). The way we compute g n (w, y, t) and 4>„(w, y, t), 
they decompose naturally as sums over the parameter (3. We define u n to be the number of 
such terms in g n (w, y, t) and v n to be the number of terms in </>^(u>, y, t). 

An examination of our construction shows that u n and v n satisfy the recursive estimates 

3 3 n n 

«n+l < Y djUn-j + Y b i V n-j + Y C l 7l "ti-j + Y C lli V n-j + (6-9) 
j=0 j=0 j=0 j=0 

3 3 n n 

V n +1 < d j U n~j + 6 J V n-j + C 3l3 U n-j + C ill V n-j, (6-10) 

3=0 j=0 j=0 3=0 

where a^, 6j, q, di, Ci and ji are fixed numbers. The exponentials 7/ arise from an estimate 
(proven in the proof of Lemma 5.2 of [15]]) for the number of Taylor series terms of any given 
order in the expansion of E(a(t) + ey). 

We substitute ( |6 . 1 U| ) for the last term in ( |6.9| ) and add the result to ( p.l(J| ). By some 



simple estimates this leads to a recursive estimate for the single quantity z n = u n + v n of the 
form 

3 n 

Z n +1 < Yl ^3 Z n-j + Y Cl 3 Zn-j- 
3=0 3=0 

An easy induction on n shows that this implies that z n grows at most like e kn for a 
sufficiently large value of k. The quantity z n is the number of terms in <f) n (w,y,t), so this 
proves the assertion. I 



Proposition now follows easily. ■ 

7 Exponential Error Bounds 



In this section, we prove Theorem 4.1. 



7.1 The Explicit Error Term 

We use the following abstract lemma, whose proof is an easy application of Duhamel's 



formula (see e.g. p~3[| ) 



Lemma 7.1 Suppose H{Ti) is a family of self-adjoint operators for h > 0. Suppose i/j(t, Ti) 
belongs to the domain of H(h), is continuously differentiable in t, and approximately solves 
the Schrddinger equation in the sense that 

where £(t, ft) satisfies 

||£(t, h) (I < n(t, h). 
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Then, for t > 0, 

|| e -ttff(a)/R^/ h s _ M t %)\\ < h- 1 f fji(s,K)ds. 

Jo 

The analogous statement holds for t < 0. 

We substitute our approximate solution ( f4.1|) 
F e iS/e2 e irra/e ( Y, e" K + e N+1 <^ +1 + t N+2 0^ +2 \ into the Schrddinger equation and 

V n=0 ) 

compute the residual term £ N . 

It is more convenient to write this term in the multiple scales notation. We also use the 
jjj(rri) (Q(f)\ ( jyi EMa(t\) 

notation e m ; y m to denote the Taylor series term e' J ' - yK 

ml , jl 

\j\=m J 

In this notation, the residual ^(w,y,t) is given, up to a phase factor, by two sums of 
terms. The first one contains all terms that do not involve derivatives of the cut-off. The 
second contains all terms that do involve derivatives of the cut-off. 
The first sum is F(w) times the following: 

e N+3 

— (A^v-i) $ (7.1) 

+ — (A w g N ) $ (7.2) 

+ e^ 3 (V w g N -i) ■ (V w $) (7.3) 

+ e N+i (V w g N ) ■ (V„$) (7.4) 

+ —9n-i (A^$) (7.5) 

e N+i 

+ — g N (A ffi $) (7.6) 

,AT+3 

+ — (A^J (7.7) 

,iV+4 

+ — (A,,,^) (7-8) 

+ e N+3 (V w ■ V y g N ) $ (7.9) 

+ e^ 3 (V y g N ) ■ (V tt $) (7.10) 

+ e^ 3 (V w -V w ^) (7.11) 

+ i^ASr+i (7-12) 

+ it N+A ^N + 2 (7-13) 
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e iV+5 



+ — (a w ^ +1 ) 



+ e N+4 (V„ • V^ +1 
+ e^ 5 (V. • V^ +2 ) 

+ — K^+i) 

e Af+4 

+ — ( A V^+ 2 ) 

h3 

■^(a(t))y 2 ^ +1 



2 

4 

E( 2 )(a(t))y 2 0^ +2 



2 

- J2e N -(E(a(t) + ey)- £ e " g^fll y" ) gJV _ n$ 

n=0 \ m<2+n m - 

~ £ e N ~ n (E(a(t) + ey) - £ e" )} y™ [ 

n=0 V m<2+n Ul - 

- [E(a(t) + ey) - £ y™ ] ^ N+1 

m<2 m - 

e N ^ \E(a(t) + ey) - £ e ^^fV ] 0^ +2 - 

m<2 m - 

The second sum arises from terms in which the cut-off F(w) is differentiated. It is 

N n+A 

£ — (A W F) g n $ 

n=0 Z 
W+2 e «+4 



+ £ V ( A - F ) *» 



n=0 



iV 



+ £ e"+ 4 (V w F)-(V^ n )$ 



n=0 



+ £ e n+4 g n (V w F)-(V w $) 

n=0 
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N+2 

+ E ^ (V«F) • (V w <^) (7.30) 

n=0 



TV 



+ E e n+3 (V W F) ■ (V^ n ) $ (7.31) 



n=0 



+ E ^" +3 (V«F) ■ (V„^) (7.32) 



n=0 



7.2 Optimal Truncation 



Each error term in the first sum ( |7.1| )-( |7T25| ) can be written as a uniformly bounded function 



times one of the following two forms: 

A = $?(w,t) E E c r,j(w,t)ifj(y,t) 

r \j\<P(r) 

& = E E d r>,j{™,t) <Pj(V,t), 
r> \j\<p'(r') 

where ^(w,t) G TC e i, VjiVit) — e ~ d ^ 2( Pj{A{t)i B(t), 1, 0, 0, y), r, r' denote a collective set of 
indices that belong to some finite set, and p(r) and p'(r') limit the number of multi-indices 
j allowed in the second sum. 

The error term £(u>, y, t) G H. e \ needs to be estimated for t G M, in the following norm 

IKWII = {/ Rd ll£ (*-«(*), (*"*(*))/*,*) \\n el dx" 1/2 

~) 1/2 

||^(w,^/e,t)||^ el ^ 



With that norm, using the Cauchy-Schwarz inequality and the L 2 (IR d ) orthonormality 
of the (pj(y,t), we obtain the following estimate for the norm of A in terms of the norm of 
vector c r (w,t) G Z 2 (N d ,(Z7), 



1/2 



< e su p ll*(<M)iki sup \\c r ( w ,t)\\ E 1 • ( 7 - 33 ) 



w&uppF uiGsuppF 



\j\<p(r) 



By similar arguments we get the following estimate for the norm of B in terms of the norm 
of the vector d r >(w,t) G / 2 (N d ,H el ), 

||B|| < E SU P \\d*(w,t)\\ ( E A ■ ( 7 - 34 ) 

r , «,esu PP F \ |i|<p'(r') / 

Note also that 

£ 1 < C' (r ] + d ). (7-35) 

\j\<p>(r>) V ' 7 

which grows at most polynomially with p'(r'). 
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Lemma 7.2 For t G [0, T], and for any a G N d and 7 G N d ; i/iere exzsi Co > and r > 0, 
snc/i that 



E E E E 

/3eSn,l P<n fc+|Z|<p+f \j\<J+n+2(p-\l\-k) 



*(tw, t)D*D^c n)P)l;k>P)j (w, t)<fj(y, t) 



< 



Co {nVVo}' 



(7.36) 



and 



E E E E 

£eB„, 2 p<n-l fc+ | Z |<p + IlzJ. y|<j +n _i + 2(p_|J|_fc) 



< C {n^ro} n 

If the operator is replaced by the operator y 1 , the same bounds are valid. 



(7.37) 



Proof: We begin with ( |7.36|) . We have 

E D w c n >P> l,k,l3,A w ^) D yfj(y,t) 

\j\<J+n+2(p-\l\-k) 



E 

|fe|<J+|7l+"+2(p-|i|-fc) 



(if, Dfo) D^Cn tPtljktP (w,t)) f k (y,t). 



We know that the vector (ip, D^ip) D™c n>Pi i t k, p(w, t) satisfies the estimate 



(tp,Dltp)D*c niPiltkjf3 (w,t) < Di D\ 



|a|+|i|+4n (a + 1)1 \t\ P k k 



X 



(l + H)^ 1 p\ 5 k 
1 (J + \-?\ + n + 2(p - \l\ - k))\ 



(\\B\\do) hl 



(7.38) 



Here 5 > is the distance in the complex plane from [0, T] to the complement of Q. 

Since the number of indices in B n 1 is bounded by e K ° n , D 2 > 1, and (a+l)\ < {\a\ + \ 
we can estimate the sum (|7.36|) by 



^(l + H)^ 1 



o K "" /;.]'" 2 

p<n 



E 



\l\+k<p+n/2 



x (\a\ + \l\)\ ^(J+\ 1 \ + n + 2(p-\l\-k))l 



(7.39) 



Then, using a\b\ < (a + b)\, the fact that (a + 2p)\/{p\) 2 is increasing in p, and p < n, we 
have 

(J+|7|+n + 2(p-|Z|-fc))! ((|a| + |/|)!) 2 ^ ( J + | 7 | + 2|a| + 3n - 2k)\ 



< 
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so that (|7.39| ) is bounded by 



grfPHg e K n D Un/2 ^ , ,„ ^ AY 

n!>/j!(H-H)^ 2 kn h W 



(J+| 7 | + 2|a| + 3n-2A;)! £ 1. 

|£|<p+n/2-fc 

The last term is bounded by ^ P 77 "/^ + ^ < aoe 3o-n/2^ w h ere | x j denotes the integer part 
of X. 



Using fc 2fc < (2k) 2k , a a b b < (a + b) a+b and a! < a a we have 

(J + |7| + 2|a| + 3n-2A;)! A; 2fc < ( J + | 7 | + 2\a\ + 3n) 
Since we can assume without loss that 5 < 1, this implies 



J+|7|+2|a|+3n 



E 

k=0 



6 



(J+| 7 | + 2|a| + 3n-2A;)! 



p+n/2 

< (J + | 7 |+2|a| + 3n) " 2 M £ <T fc 

k=0 

„ f T II r.1 I o , J+l7l+2|c|+3n Jf-j 

< (J+| 7 |+2|a| + 3n) a ^ 5 p , 

for some constant A'i that satisfies <5 _1 — 1 > iff" Together with 



(7.40) 



E (t/tr < 

p<n 



K 2 (t/5) n if t/S > 1 
fsf™ if t/S = 1 
£T 2 if t/S < 1, 



(7.41) 



where is constant, we get (in the first case above) 

E E E I D ]f) DZ c n, P ,i,kA w ^) 

PeBn,! p<n k+\l\<p+% 

a K 1 K 2 D 1 D 2 a \\\B\\d )^ e ^ +3CT / 2 )"Z) 2 lln/ V 



< 



J+|7|+2|q|+3n 



(J+ | 7 | +2|a| +3n) : 5 



(1 + H) rf+1 y/j\n\S 3n / 2 

We postpone the study of the dependence of our estimates on t and J to Section ||. So, using 
the above, 

J+|7|+2|a|+3n J+|-v|+2|a| 3n 

(J + | 7 | + 2\a\ + 3n) < (J + | 7 | + 2|a| + 3n)^^((J + | 7 | + 2\a\ + 3)n)~ , 

and the existence of < a < b, such that a n n n < n! < 6 n ra n , we learn the existence of 
positive constants (i.e., independent of n) K 3 , K4 and -fr 5 , such that 

3n/2 



E E E || <^£>> || < K ^ K t^ < K 3 (K 5 n^r. 

peB n> i P <n k+\i\< P +% a n 
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This yields the result with Co = K% and tq = K§. 

The second sum is dealt with in the same manner, since the vectors d n ^i^,p{w, t) satisfy 
the same bounds as c ntPt i t k,p(w,t) does with n replaced by n — 1. 

Finally, the replacement of by the operator y 1 means that the matrix (ip\D^ip) must 
be replaced by the matrix ((ply^ip). But the latter has the same properties as the former; 
the bounds above remain true with replaced by \\A\\ from ( [7.38 ) onward. This affects 
the definition of C Q only. I 

The following lemma is the key to the proof of exponential accuracy of our approximation 
by means of optimal truncation. 

Lemma 7.3 For sufficiently small g > 0, there exist T(g) and C(g) > such that the choice 
N(e) = [gVe 2 ] implies that the norm of the error term £jv( e )(£) given by ^T\ ) satisfies 

UN(e)(t)\\ < C(g)e-^\ 

Proof: The previous lemma, formulas Q7.33 ), ( 7.34 ) and ( |7.35|) show that all terms in the 
first sum defining £jy except ( |7.12|) , (|7.13|) , (|7.22|) , and (|7.23| ) are exponentially small, once 
we prove 

C e N U (iV(e) 1 / 2 r } m < Ce- r ^. (7.42) 

Because g 2 /e 2 — 1 < iV < g 2 /e 2 , if we choose < g < l/r , the left hand side of this 
inequality is bounded by 

C {eN 1 / 2 ^ < C {gr } N < C e" 1 ln(sTo)|Af < C e |ln(9ro)l e ~l i-(^o)| 5 2 A 2 ? (7.43) 

which gives 

C(g) = C e\ 1 ^ and T(g) = |ln(^r )U 2 . 
The terms ( |7. 12| ) and ( [7. 13| ) can be dealt with in a similar fashion once we have computed 

E E E E d N+hP:lAP j(w,t)ip j (y,t) 

/3eBjV+i l2 P<N k+ \i\< p+ N \j\<J+N+2(p-\l\-k) 

+ d N+ltPj i tk:l3 j(w,t)'f j (y,t), 

where the second term equals 

E E E E ((U^\v>) 

/3£B N+1>2 p<N fc+ |;|< p+ A[ \k\<J+N+2+2(p-\l\-k) 

iEW(a(t)) , 2 ,\ . . \ 

Lemma [67^ shows that £Z/v+i )P) j 1 fc )( 3 satisfies bounds similar to those satisfied by d,N+i, P ,i,k,p,3 
and the term above is taken care of by lemma |7T2j Similar statements are true for <fijj +2 , and 
the analysis above also applies to these error terms. 
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Next consider ( [7.221 ). By the mean value theorem, there exists ( q (y,t,e) = a(t) + 
q (y, t, e) ey, where q G N d and q (y, t, e) G (0, 1), such that 



E(a(t) + ey)- £ e m g^Q) y 



m<2+n 

Hence, we need to estimate 



ml 



E 



D*E(( q (y,t,e)) 



|<?|=2+n+l 



N 



E E 



D«E(( q {y,t,e)) 



(7.44) 



n=0 |q|=2+n+l 

with the following restrictions 

\j\ < J+(N-n) + 2(p-k-\l\) 
k+\l\ < p+(N-n)/2 
p < N — m 

[3 G Bi^N-n- 

We take a fixed value of n G [0, N], and consider the vectors 

D«E(( q (y,t,e)) , 

j {if, y q (p){t)c N - ntPt i tkt0 {w,t) 



we have to estimate. Due to the presence of the cut-off function F (which we have omitted 
in the notation), we have 



(7.45) 



\DiE(( q (y,t,e))\ < c 6 D 2 



\q\ 



(i + M) 



<f+l : 



and with our bounds on the matrix (ip, y q ip)(t) and on the vector cjj-n,p,i,k,p{w,t), we can 
write 



D«E(( g (y,t,e)) 



(7.46) 



< 



c 6 Z^+ 3 (d \\A\\) n+3 J(J + 3 + N + 2(p-k 



D l D 



l\+4(N-n) 



\f\P Uk 



(1 + (n + 3)) rf+1 v/j[ 1 2 ' p! 5 k * 

Then we use similar estimates to the above and the restrictions ( 7.45 ) to get 



k 2k l\l\ 



(J + 3 + N + 2(p-k- 



p\ pi 



- < (2k) 



2k (J + 3 + N + 2(p~k))\ 



< (2k) 



2k 



(J + 3 + 3iV - 2n - 2k)\ (J + 3 + 3N-2n) 



pi pi 

J+3+3N-2n 



< 



(N-n)\(N-n)\ ~ (N — n)\ (N — n)\ 

Using this and \l\ < 3(N — n)/2, we see that (|7.46|) is bounded above by 

n llJV/2 



4 d + 1 VJ! 



jfn/2 (jv-n)!<j* 
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(j + 3 + 3Ar_2n)( J+3+3iV )/ 2 . 



Finally, with 



N-n PH — 

E E 

p=0 fc=0 



N-n 



(see O), (CT ), the bounds E\i\<p+(N-n)/2 1 < a e 3 ^"")/ 2 , E\ q \< n+3 1 < 
and |jBi,jv-n| < e - K 'o(w-n) ) we g e ^. (.^h conditions ( |7.45| ) on the summations) 



<j(n+3) 



E E < 

|g|=2+n+l P,p,k,l,j 



N 



D«E(( q (y,t,e)) 



(<p, y q ^){t) c N _ ntP> i >k;/3 (w,t) 



< 



^ K i K 2*DjPMM) 3 e5 , N /2 d n ( do \\ A \\f e N ( J + 3 + 3iV)( J + 3 + w )/ 2 



D 



9/2 



7V-n 



X 



5 1/2 doll^l 



*\^~ n (7 + 3 + 3^)^-" 
5/ (N-n)\ ' 



(7.47) 



Postponing the study of the t and J dependence of our estimates, we use the bound 

( J + 3 + 3iV)( J + 3 + JV )/ 2 < N N ' 2 ( J + 3 + 3iV) (J+3)/2 (J + Q) N/2 to establish the existence of 

constants L ,Li,L 2 , independent of N and n, such that (|7.47|) is bounded above by 



L Q L^e N N N ' 2 



(L 2 N) 



N-n 



(N-n)\ ' 

It remains for us to sum over n and use ( |7.33|) to bound (|7.44|) by 

[L 2 N) N ~ n 



N 



e 3 v ^- e 3 ^ 2 L L^e N N N / 2 £ 

n=0 

(X 



(L 2 N) S 

s=0 

< e 3 v ^- L (e^ 2 L 1 e L2 fe N N N / 2 . 



si 



If we choose g < l/(Li e i2+3ff/2 ), we can apply the analysis ( [7.431 ) to obtain an exponentially 
small bound on ( [7.441 ) by the optimal truncation N(e) = [g 2 /e 2 ]. 

Since the estimates we have on the ef s are similar to those we have on the c's, with the 
replacement of n by n — 1, the same exponential bound is valid for (|7.23|) , (see (|7.34|) ) and 
the analysis of the the first collection of error terms is completed. 

We now need to take into account the error terms fl7jp to ( 17732]) arising; from the 
derivatives of the cut-off function F. Choose Fq > that satisfies 

max {\A w F(w)\, \\V w F(w)\\} < F , 

uniformly in w, and recall that for any % = 1, . . . , d, 

supp d Wi F(w) C {w eM d : b < \w\ < b x ] (7.48) 
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for some < 60 < b\ < 00. Now consider ( |7.26| ). We express g n in terms of the c's to see 
2 

that the norm of - times ( ff^BQ (in L 2 (R d ,H c i)) can be bounded as follows: 



N 



A w F(w)g n (w,y,t)e n $(w,t) 

n=0 



N 



< E e n x , 



A w F(w) c nJ (w,t)e- d / 2 (p j (w/e,t) 

\j\<J+3n 



\Mw,t)\\ 2 Hei dw 



N 



< F J2t n sup \\c n (w,t)\\ 

n=0 tuesuppFCR 1 * 



J2 / \e- d / 2 <pj(w/e,t)r dw. (7.49) 

\ \j\<J+3n J \ w ^ bo 



We know from Section 7 of |15[] that there exists a constant < /3d depending on the 
dimension d only, such that 



\/m + d < W(N|e), for all |j| < J + 3iV 
and I j I < J + 3 iV imply 



|tu|>bo 



-*/ 2 <Pj(w/e,t)\ 2 dw < e &b1 e -(^)/(i2||A|P^)_ 



All the conditions here will be satisfied if N(e) = [g 2 /e 2 ], provided we choose g and e to 
satisfy 

e 2 (d + 2J) + 6g 2 <b 2 /\\A\\ 2 . 
For such a choice, using Y,\j\<j+&n e 2/3d l J l < cr e^ +2/3d ^ J + 3 ™) ; we get 



£ / |e-<VV>A,0| 2 ^ < ^ e ( CT+ 2 /3d )(J + 3.)/2 e _ (6g)(12 ||AP^)_ 

^ |j|<J+3n J M>bo 

Moreover, by means of manipulations that by now are familiar, 

\\Cn(w,t)\\ < Yl E E 11^,^/3(^^)11 
/36B„,i P<n fc+|«|<p+f 

l\+4n„ \t\ p k k (J + n + 2(p-\l\-k)) 



< E E E d 1 d^i\ 



PeB n ,i p<n k+\l\<p+% 



< E E E 

/3eB„,i p<n fc+|Z|< p+ 2 



p! <5 fc 



V! 



D 1 D 1 2 ln/2 \t\P (J + 3nY J+ 



3n)/2 



n! 



< e*»» a e 3 ™/ 2 J^Jfc 1 (M) " ( j + 3 n)^ 2 



(J + 3iV) 3n/2 



n! 
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Combining these estimates, we get the existence of positive constants M and Mi, such that 
foriV = [g 2 /e% 



N 



& w F(w) g n {w, y, t) e n t) 



n=0 



< e -(6g)/(12||A||^) J2 M i eM l N3/2 T 
n=0 



III 



< e -(62)/(12||A|] 2 e 2 ) MQe 6MiJV 3 /2 < e -(^)/(12||A|| 2 e 2 ) Moe M lS 3 A 2 < Mq e ~ (ft§)/(24||A|| V) ? 

provided 

M^ 3 < 6g/(24||A|| 2 ). 

All other terms in the list (fT26|) to (fT3^) can be estimated in a similar fashion under a 
similar condition on g. 

This concludes the proof of our lemma. M 

Remark: It is not difficult to check that if we keep N fixed, then our approximation ( |4.1| ) 
ip(w,y,t) is accurate up to an error of order e N , as expected. 

A by-product of our estimates on the terms stemming from the introduction of the cutoff 
is that our approximation is exponentially localized in a ball centered at a(t) of any radius 
60, as stated in the second part of Theorem |4.1|. 



Hence, we have completed the proof of Theorem O 



8 Generalizations 



As in ||15|| , under some mild supplementary assumptions, we can extend our results to allow 
< t < T(e) with T(e) ~ ln(l/e 2 ). This proves the validity of our construction up to the 
Ehrenfest time scale. 



Theorem 8.1 In addition to the assumptions of Theorem \4-l\ , assume that a classical solu- 
tion to the equation ( $-4 ) exists for all t e M. Moreover, assume that for all z in a complex 
neighborhood of S ; the following bound is satisfied 

\E(z)\ < N e M ' 2 ', 

and that E(x) is bounded below. Suppose also that there exist L and A > ; such that for all 
t e R 

\\A(t)\\ + \\B(t)\\ < L e xt . 

Then, there exist r 1 , C , T' > 0, and < a, a' < 2 such that the approximation defined 
by choosing N(e) ~ l/e a is accurate up to an error whose norm is bounded by C ' e~ T ' ^ a , 
uniformly for all times <t <T' ln(l/e 2 ). 

Proof: It is enough to mimick the proof of the corresponding result for the semiclassical 
propagation of the Schrodinger equation in ||15|| , since our hypotheses imply that nothing can 



happen on the adiabatic side of the problem. By the conservation of energy, the exponential 
bound on E(z) and the assumed existence of a Liapunov exponent, we easily see from the 
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proof of Lemmas |7.2| and [7.3| , that the behavior in t of all constants (independent of N) is 
at worst exponential in t. From the conditions D 2 > e KT , with K some constant, we need 
to take g(T) < goe~ 9lt so that the optimal truncation procedure yields an error of the order 
e K T e -g*e-^i/e^ The choice T ^ < T'ln(l/e 2 ), with V > sufficiently small, gives the 
desired result. I 

Similarly, we can extend our results to allow initial conditions in a wider class of vectors. 
Indeed, we have been careful to make explicit the J dependence in all estimates so that 
we can control the error term as a function of J. Recall that J is fixed arbitrarily in ( |3.8| ) 
which gives the expansion in the basis <pj(A(Q), B(0), e 2 , a(0), T](0),x) of the nuclear part of 
the wave function that we take as an initial condition. 

As in [HJ, for (a, 77) G IR 2d , we introduce the operator A e (a, rj) such that 



(A £ (a, V )f)(x) = c-*e^ f((x-a)/e). 
We define a dense set C in L 2 (IR d ), that is contained in the set S of Schwartz functions, by 

C = < /(a;) = E c 3 L Pj$^ 1) 0; 0; x ) e <5> sucn that 
I J 

there exists K > with V] |cj| 2 < e~ KJ , for large J 

lil>J 



Remark It is easy to check that the inequality in ( |8.1| ) is equivalent to the requirement 
that the coefficients of / satisfy 

\ Cj \ < e- K W 

for large Another equivalent definition of C is 

C = U t>0 e- tHh ° S, 

where Hho = — A/2 + x 2 /2is the harmonic oscillator Hamiltonian. The set C is also called 
the set of analytic vectors pTj for the harmonic oscillator Hamiltonian. 



Let / e C. We set 

fj(y,t) = E c lVj (A(t),B(t),e 2 ,0,0,y), and 
\j\<J 

j 

where the classical quantities a(t), rj(t), A(t), B(t), and S(t) correspond to the initial con- 
ditions a(0), 77(0), A(0) = B(0) =E, and 5(0). We consider the construction described in 
Section [| corresponding to the initial condition g o (0,y,t) = fj(y,t), making explicit the 
dependence on J in the notation: 

^j,N(w,y,t) 

/ N N+2 \ 

= F(w) e iS ^ 2 eWW' E ^ 9n,j(v>, V, t) $(w, t) + E ^ <jK V,t) )■ 

\n=0 n=2 / 
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Recall that 



V JtN {w,y, 0) 

/ N+2 \ 

= F(w) e iS ^ 2 /j(y, 0) 0) + £ e"<^>,y,0) . 

V n=2 / 

Let 1/ > 0, and consider N(e) = {g 2 /^ 2 } and J(e) = uN(e). We define our more general 
initial conditions as 

*f{w,y, o) 

/ AT(e)+2 

= F(w) e iS(0)/e2 e ir,{oyy/e f(y, 0) 0) + £ e n < J(e) K 2/, 0) 

\ n=2 

which corresponds, when we get back to the variables (X,t), to an initial state *£?f(X — 
a(0), (X — a(0))/e, 0) whose projection along the electronic eigenvector 0) yields a 

nuclear wave packet of the form (A e (a(0), T](0))f)(X). Note that the component of the 
initial state perpendicular to $(X, 0) necessary to achieve exponential accuracy depends on 
e. This component is determined by the coefficients of the function /. 
We can now state our result for such general initial conditions 



Theorem 8.2 Assume the hypotheses of Theorem ^ ■ 1\ and consider the above constructions. 
There exist sufficiently small g > and positive constants C(g), T(g), such that with the 
definition 

t, e) = $ Ae)Me) (X - a(t), (X - a(t))/e, t), 

we have 

\ e -*HW*Hf f ( X ,0,e) - MX,t,e)\\ LHRdtHei) < C{g)e-^\ 
for allt e [0,T], ase^O. 

Moreover, the result for times T ~ ln(l/e 2 ) corresponding to Theorem \8.1\ is also true for 
these initial conditions. 

Proof: We have 

e -itH(e)/e* % f (X,0,e) 
= e -itH(e)/e* 0> e) _ 0> e)) + e -iW(e)/e* ^ ^ 

= **(X,t,e) + O(||e-^ 2 vI/,(X,0,6)-vI/,(X,t,6)|| L2(]RdiWci) ) 

+ O ( \\* f (X, 0, e) - 0, e) || L2( ^, Wel) ) • 

By our choice of function /, the last term is exponentially small in 1/e 2 . The remaining 
norm to estimate corresponds to the situation of Theorem ^J] in which we let the parameter 
J grow as 1/e 2 , according to our choice of J(e). But, as in the proof of Theorem 3.6 in |L5| 
for the corresponding result in semiclassical dynamics, we have made the dependence in J of 
all the key estimates explicit. It is enough to go through the proof of theorem JO] to check 
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that with J = uN, all arguments can be repeated to get the same N and e behavior for the 
estimates on the error terms, (see |15| for details). Hence, we see that for sufficiently small 
g, we can approximate the solution corresponding to these generalized initial conditions up 
to an error of order e~ r (s)/ e The Ehrenfest time regime is dealt with similarly. I 



9 Technicalities 

In this section we give the proofs of the auxiliary lemmas we used in the course of the main 
argument. 



Proof of Lemma |6.1| : We first consider the case fc > 1. By Cauchy's formula, we can 
write 



9{s 



2ni Jr (t — s)' 



ds, 



(9.2) 



1 



where T is the circular contour with center t and radius (5 — |Imt|). 



For s on T, we have (5 — |Im s\) > 



fc + 1 

< C k k (5 - \\ms\)- k < C k k 



fc + 1 
(5 — \Imt\). Thus, 

k 



fc + 1 



(S - \lmt\) 



So, by putting the norm inside the integral in (|9.2|), we have 



k 



fc + 1 

C (fc + l) fc+1 (5 - Ibntl)-*- 1 . 



{6 - \lmt\) 



fc+1 



(5 - \lmt\) 



-2 



For fc = we use the same argument with the radius of F replaced by a (5 — |Imt|) for 
any a < 1. This yields the bound 

\\g'{t)\\ < Ca- 1 ^ - Ibntl) -1 . 

The lemma follows because a < 1 is arbitrary. H 



Proof of Lemma |6.4| : To prove the quantity v is finite, we estimate 

1 1 



E 



{l :0< k < ai} + (l+l«-/|) d+1 

V - 1 1 

^ n -i- 

{ l : < U < oii } 

H 

2 



<[¥] 



(l + |/|) d+1 (l + |a-/|) d+1 
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+ E 

{ / : < k < a, } 



(l + \l\) d + l {l + \a-l\) d + 1 



< 



< 



2<i+2 



1 + Id 



\d+l 



E 

{ l : < Z< < a, } 



\d+l 



2 d+2 



1 + a 



\d+l 



(l + l'l) 



Thus, i/ < 2 d+2 £ (1 + IZI)-**- 1 . 



To see that the right hand side of this inequality is finite, we note that the number of multi 
indices / with 1/1 



L is the binomial coefficient ( ^ i ^ ) ' w ^h the convention that 








1. Thus, 



v < 2 d+2 Y, 



L=0 



L + d-1 
d-1 



0- + L) 



d+l 



2^+2 



(L + d-l)(L + d-2)---(L + l) 



For large L, 



(L + d-l)(L + d-2)---(L + l) 



(L + l) 



d+l 



is asymptotic to L , so v is finite. 



Since D a (M N) = V 

{l:0<k< ai } 



[] ^ " j ^ (D'M) (D^-')iV), we have 



< 



£> a (MiV))(x) 



s n n 

{l:Q<k<ai} \_j=l \ 3 



mix) n(x) a(x)\ a+ ^ (/ + p)! (a - I + q)\ 



= m(x) n(x) a(x) ]a+p+ql (a + p + q)\ 
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x E 

{l:0<k<a t } 



n 



Q i 



3=1 \ 



atj + pj + qj 

h + Pi 



Since + * 

V h + Pi 



~ { I, 



(l + \l\) d + 1 (1 + \a-l\) d + v 
we therefore have 



(D a (MN))(x) || < m(x) a(x) |a+p+g| (a + p + g)! 



1 



X , (l + |/|) d+1 (l + |a-/|) rf + r 



< m{x) n(x) v a(x) 



{l:0<h< ai } 

\a +P+q \ (a+p + g)! 



;i + |a|) d + 1 ' 



Proof of Lemma |0|: If f(t) satisfies < C \t\ p dist(t)" fc , for all t e £1, there exists 

g(t) analytic in Q, such that f(t) = t v g(t) and ||g(t)|| < G dist(t)~ fc . We use the integration 
path from to t G Q parametrized by j(u) = tu, with u G [0, 1], to compute 



f(s) ds 



f(tu) du 



t (tu) p g(tu) du 



< C\t 



P+l 



it?- 



du 



< ( ' i-i-j dist(t)- 



/o dist(tM) fc p + 

since, by assumption, dist(wt) is a decreasing function of u 



(9.3) 
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